www.gusucode.com > matlab编程NSCT分解 图像融合 各个融合指标评价体系 分解源码程序 > matlab编程NSCT分解 图像融合 各个融合指标评价体系 分解源码程序/NSCT/dfbdecdemo.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Image decomposition by nonsubsampled contourlet transform (NSSC). % This is the iterated filter bank that computes the nonsubsampled % contourlet transform. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Parameteters: nlevels = 4 ; % Decomposition level %pfilter = 'pkva' ; % Pyramidal filter dfilter = 'dmaxflat7'; %'cd' ; % Directional filter % Nonsubsampled Contourlet decomposition im=imread('C:\Users\Administrator\Desktop\Risa-f4c8a9f7c3a06c83224db118c2ef871234d3f6e6\nsct_toolbox\zoneplate.png'); coeffs = nsdfbdec( double(im), dfilter, nlevels ); disp( nlevels); disp(dfilter); % Display the coefficients %disp('Displaying the contourlet coefficients...') ; %shownsct( coeffs ) ; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Nonsubsampled Contourlet transform (NSSC) reconstruction. % This is the inverse of nsscdec, i.e. % imrec = nsscrec(coeffs, dfilter, pfilter); % would reconstruct imrec = im %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Reconstruct image imrec = nsdfbrec( coeffs, dfilter ) ; % Show the reconstruction image and the original image %if 0 figure; subplot(1,2,1), imagesc( im, [0, 255] ); title('Original image' ) ; colormap(gray); axis image off; subplot(1,2,2), imagesc( imrec, [0, 255] ); title('Reconstructed image' ) ; colormap(gray); axis image off; %end mse = sum( sum( (imrec - double(im)).^2 ) ); mse = mse / prod(size(im)); disp( sprintf('The mean square error is: %f', mse ) ); disp(' ');